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YEAR  1 : 

Task  1 :  Development  of  experimental  approaches  and  their  adaptation  for 
gas/surface  systems  declared  in  the  project. 

A  hot  filament  method  is  used  to  study  the  heat  transfer  between  tungsten  surface  and 
hollow  cathode  glow  discharge  plasmas  in  argon  and  C02.  The  dependence  of  the 
electric  power  supplied  to  a  tungsten  wire  on  the  discharge  current  is  determined  for 
argon  and  carbon  dioxide  in  the  temperature  range  between  1000K  and  1700  K.  A 
difference  in  heat  transfer  at  the  tungsten  wire  surface  is  found  between  experiments  on 
argon  and  carbon  dioxide.  The  difference  is  attributed  to  heterogeneous  recombination 
in  C02  plasma. 

It  is  hardly  possible  to  enumerate  all  applied  and  fundamental  issues  related  to  energy 
transfer  across  plasma-solid  interfaces.  Problems  of  this  kind  arise  in  areas  ranging 
from  aerothermodynamics  to  plasma  enhanced  chemical  vapor  deposition,  sputtering, 
and  etching  technologies.  Therefore,  research  on  this  subject  is  important  from  both 
theoretical  and  practical  perspectives.  Of  particular  importance  are  further  experimental 
studies  using  various  techniques  with  a  view  to  improving  theoretical  models  of  plasma- 
surface  interaction. 

The  choice  of  a  hot  filament  technique,  widely  used  in  studies  of  heat  transfer  in  low 
pressure  gas-surface  systems  starting  from  Knudsen’s  work,  has  the  technical 
advantage  of  simple  processing  control. 

Tungsten  is  used  as  a  filament  material  for  the  following  reasons.  First,  tungsten  is  a 
metal  with  well-known  electrical  and  thermal  properties.  Second,  it  has  a  high  melting 
point  of  3653  K,  which  makes  it  possible  to  conduct  investigations  in  a  wide  temperature 
range.  Moreover,  a  vast  amount  of  data  is  available  on  heat  and  mass  transfer  at 
tungsten  surface. 

Plasma  is  produced  in  a  hollow  cathode  glow  discharge.  A  key  advantage  of  this  system 
compared  to  those  with  planar  electrodes  is  the  feasibility  of  a  high  current  density 
through  the  electrode  gap.  Furthermore,  it  can  be  used  to  produce  highly  ionized 
plasmas  at  low  gas  pressures  and  provides  good  spatial  resolution  of  glow  discharge 
zones. 

Argon  and  Carbon  dioxide  are  used  as  working  gases  because  they  have  similar 
molecular  masses,  but  widely  different  physical  and  chemical  properties  determining  the 
characteristics  of  processes  that  occur  both  in  the  gas  phase  and  at  gas-solid 
interfaces. 

Figure  1 .1  shows  a  block  diagram  of  the  experimental  setup.  A  glow  discharge  chamber 
is  mounted  on  insulators  inside  a  vacuum  chamber.  The  vacuum  chamber  is  evacuated 
to  a  residual  gas  pressure  on  the  order  of  10“3  Pa  by  oil  diffusion  pump.  The  vacuum 
chamber  is  a  410  x  430  x  500  mm3  parallelepiped  made  of  the  stainless  steel 
12X18H10T. 

The  working  gas  is  fed  into  the  discharge  chamber  through  an  inlet  duct  of  3  mm  in 
diameter  via  mass  flow  controller  and  evacuated  through  a  9  mm  diameter  outlet  in  the 
opposite  flange.  A  gas  flow  rate  set  at  approximately  0.2  cm3  atm/s  maintains  a 
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discharge  chamber  pressure  of  about  1  Pa  when  the  vacuum  chamber  pressure  is 
approximately  4  x  10“2  Pa. 


Fig.  1.1.  Experimental  setup:  (1 )  vacuum  chamber;  (2)  gas 
discharge  chamber;  (3)  tungsten  wire;  (4)  gas  inlet  system; 

(5)  filament  heating  and  temperature  stabilization  system; 

(6)  evacuation  system;  (7)  insulators;  (8)  discharge  power 

supply. 

Task  2:  Surface  preparation  and  diagnostics  by  AFM. 

The  surface  preparation  procedure  depends  on  the  task.  The  same  is  for  the  surface 
diagnostics.  In  particular,  in  experiments  devoted  to  surface  composition  effect  on  gas 
flow  rate  in  channels  the  system  for  exposure  of  evaporated  metal  or  adsorbed  gas  is 
used  (see  below  Task  7).  In  experiments  devoted  to  surface  chemical  composition 
effect  of  heat  exchange  between  rarefied  gas  and  a  solid  body  the  surface  preparation 
procedure  mainly  includes  purification  cycles  based  on  high  temperature  heating  (see 
Task  6).  In  both  cases  Auger  Electron  Spectroscopy  as  a  means  of  surface  diagnostics 
is  used. 

The  problem  of  surface  stricture  effect  on  character  on  gas/surface  interaction  and  as  a 
result  on  accommodation  coefficients  needs  attracting  such  means  of  diagnostics  as 
Atomic  Force  Microscopy  (AFM).  In  particular,  for  Si-based  samples  the  following  data 
have  been  obtained  with  the  use  of  AFM  technique  (Fig. 2.1 ,  Fig2.2  and  Fig. 2. 3). 
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Fig. 2.1 .  Surface  structure  for  Si-based  material  according  to  AFM  data:  3-dimensional 
image. 


X[Mm] 


Fig. 2. 2.  Surface  structure  for  Si-based  material  according  to  AFM  data:  separate 
profiles  along  one  scanning  direction. 
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Fig2.2.  Surface  micro  roughness  distributions  function  for  investigated  sample 
according  to  AFM. 

These  data  were  used  for  simulation  of  surface  structure  and  calculation  of  gas 
flow  rate  in  channels  (see  Task  5). 


TASK  3:  NORMAL  MOMENTUM  ACCOMMODATION  STUDY. 


A  dynamic  technique  for  normal  momentum  accommodation  study  in  a  gas-surface 
system  is  developed.  The  method  can  also  be  used  for  pressure  measurements  and 
adsorption  studies. 

Characteristics  of  dynamic  pressure  gauges  allow  their  application  in  various  fields  of 
science  and  engineering.  For  example,  by  measuring  the  frequency  shift  of  sensor 
oscillations,  one  can  measure  gas  adsorption  on  the  sensor  surface.  Another  important 
application  of  the  dynamic  technique  is  the  investigation  of  momentum-transfer 
efficiency  in  the  course  of  interaction  of  the  gas  and  solid  surface,  which  is  generally 
characterized  by  accommodation  coefficients.  The  values  of  accommodation 
coefficients,  in  particular,  the  normal  momentum  accommodation  coefficient  (NMAC), 
are  mainly  obtained  by  the  molecular  beam  technique.  Implementation  of  this  method, 
however,  faces  significant  difficulties  including  the  problem  of  interpreting  the  results 
and  using  them  to  describe  the  processes  of  gas-surface  interaction  under  weakly  non 
equilibrium  conditions. 

The  proposed  method  for  NMAC  measurement  has  significant  advantages  over  other 
methods;  one  of  advantages  is  the  simplicity  of  its  implementation.  As  for  interpretation 
of  the  data  obtained,  it  is  also  rather  simple  if  some  assumptions  are  used  in  a  particular 
experimental  situation. 

The  essence  of  the  technique  proposed  is  the  measurement  of  force  acting  on  a  plate 
moving  in  a  gas.  The  plate  motion  can  be  ensured  by  free  oscillations.  In  the  present 
work,  we  measured  the  damping  coefficient  of  free  linear  oscillations  of  an  elastic  plate 
in  the  free-molecular  regime  with  varied  gas  pressure.  The  measurement  path  of  the 
experimental  setup  is  schematically  shown  in  Fig.  3.1 .  The  sensor  is  fabricated  from  a 
plate  of  microcontoured  single-crystal  silicon  by  the  MEMS  technology.  This  allowed  an 
almost  perfect  attachment  of  the  elastic  part  of  the  plate  5x5  mm  and  obtaining  rather 
stable  damping  coefficients  of  plate  oscillations  (3  =  (0.020  ±  0.002)  sec’1  under 
conditions  of  ultrahigh  vacuum  (of  the  order  of  1 0”7  Pa).  The  thicker  part  of  the  plate, 
which  does  not  directly  participate  in  the  momentum-transfer  process,  is  fixed  by 
assembling  screws  used  in  ultrahigh  vacuum  engineering.  The  elastic  plate  together 
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with  the  fixed  electrode  forms  an  electric  capacitor  of  about  15  pF,  which  is  connected 
to  the  oscillating  LC  generator  with  a  resonance  frequency  of  about  10  MHz.  The  initial 
oscillations  of  the  elastic  plate  are  excited  by  the  electric  field  between  the  sensor  and 
fixed  electrode  by  means  of  a  specially  designed  system  of  self-induced  oscillations. 
The  excitation  system  serves  as  a  chain  of  feedback  and  a  source  of  external  energy. 
The  natural  oscillation  frequency  of  the  elastic  plate  is  about  1250  Hz. 


Fig.  3.1.  Measurement  path  of  the  experimental  setup:  1)  fixed  electrode;  2) 
microcontoured  silicon  plate;  3)  high-frequency  LC  generator;  4)  frequency  detector;  5) 
selective  amplifier;  6)  oscillation-generation  system;  7)  analog-to-digital  converter;  8) 
computer. 


Fig.  3.2.  Vacuum  chamber  and  gas-inlet  system:  1)  vacuum  chamber;  2)  monopole 
mass  spectrometer;  3)  vacuum  pressure  gauge;  4)  sputter-ion  pump;  5-9)  vacuum 
valves;  10)  backing  pump  system;  11)  gas-inlet  valve;  12)  optical  membrane  gauge;  13) 
thermocouple  vacuum  gauge;  14)  gas  reservoirs. 

The  microcontoured  silicon  plate  and  the  fixed  electrode  are  mounted  vertically  in  the 
vacuum  chamber,  which  is  shown  in  Fig.  3.2  together  with  the  gas-inlet  system.  The  gas 
is  evacuated  from  the  vacuum  chamber  by  a  sputter-ion  pump  with  a  pumping  velocity 
of  100  liters/sec.  The  residual  pressure  in  the  vacuum  chamber  is  measured  by  a 
Penning  vacuum  gauge  and  reaches  less  than  10”7  Pa.  The  gas  composition  is 
controlled  by  a  monopole  mass  spectrometer  for  pressures  lower  than  10~3  Pa.  The 
vacuum  chamber  is  connected  to  the  gas-inlet  system  through  a  valve.  Preliminary 
evacuation  is  performed  by  the  backing  pump  system  consisting  of  two  adsorption 
pumps  and  a  mechanical  vacuum  pump.  The  gas  pressure  in  the  inlet  system  is 
measured  by  an  optical  membrane  gauge  and  a  thermocouple  vacuum  gauge. 
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Measurement  of  the  damping  coefficient  of  free  linear  oscillations  of  the  elastic  plate 
with  varied  gas  pressure  reduces  to  registration  of  the  value  of  a  capacitor  formed  by 
the  elastic  plate  and  the  fixed  electrode  and  included  into  the  oscillating  circuit  of  the  LC 
generator. 

For  NMAC  determination,  the  measurements  are  performed  with  fixed  values  of  gas 
pressure.  Except  for  the  standard  procedure  of  cleansing  after  the  final  stage  of 
fabrication,  the  sensor  surface  received  no  other  treatment. 

The  prior  measurements  based  on  the  use  of  this  experimental  approach  have 
demonstrated  significant  change  in  damping  process.  In  particular,  the  frequency  shift 
for  wide  range  of  Knudsen  numbers  has  been  observed  that  confirms  high  sensitivity  of 
this  method  to  gas/surface  interaction  parameters  (Figs. 2. 3,  2.4). 


Fig. 3. 3.  Frequency  shift  as  a  result  of  gas  pressure  change  (hard  line  -  1  Pa,  dotted  line 

-lOPa). 


A,  units 


Fig. 3. 4.  Decay  curve  and  its  logarithms  at  10  Pa. 

To  take  into  account  adsorption  influence  the  basic  QCM  principles  are  used. 

According  to  Sauerbrey’s  assumptions  regarding  the  mass/frequency  shift  relation  the 
density  of  deposited  thin  layer  and  elastic  modulus  are  such  that  it  can  be  treated 
approximately  as  an  extension  of  the  silicon  crystal.  The  mass  of  the  adsorbed 
molecules  Am  would  be  equivalent  to  an  increase  in  the  mass  m  of  the  elastic  plate. 
Thus  the  resonator  thickness  11  is  connected  with  its  mass  as 

h  =  —  ,  (3. a) 

Fp 

where  F  is  an  elastic  plate  square,  P  is  a  density  of  the  silicon  single  crystal  plate.  The 
resonant  frequency  of  unloaded  square  plate  can  be  presented  as 

F  =  Kh ,  (3.2) 

where  K  js  the  constant  depending  on  the  Poisson’s  ratio,  the  Young’s  modulus,  the 
elastic  plate  square  size,  the  density  of  the  silicon  single  crystal. 
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Thus  a  frequency  change  ^  can  be  connected  with  a  thickness  change  by  the 
following  correlation 


^  =  (3.: 

/  h 

Taking  into  account  the  earlier  obtained  expression  for  the  normal  momentum 
accommodation  coefficient  a"  the  decay  time  constant  P  can  be  presented  finally  as 


P  = 


(2ra»\r  W* 

ph( l  +  f)  V  nkT  ' 


(3.4) 


where  P  is  a  gas  pressure,  Mg  is  a  gas  molecule  mass,  k  is  the  Boltzmann  constant, 
T  is  a  temperature. 


The  preliminary  measurements  of  time  frequency  dependence  and  damping  process  of 
free  elastic  plate  vibrations  have  carried  out  at  ultrahigh  vacuum  conditions.  The 
temperature  dependence  of  frequency  and  spectrum  of  vibrations  are  presented  in 
Fig. 3. 5.  The  spectrum  is  obtained  as  a  result  of  the  Fast  Fourier  Transform  (FFT) 
application  to  signal  processing  at  the  rate  of  104  samples  per  second. 
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A, 

units 


T,K 


f,  Hz 


(a)  (b) 

Fig.  3.5.  The  temperature  dependence  of  natural  frequency  (a)  and  frequency  spectrum 
(b)  of  the  elastic  plate  vibrations:  A  is  amplitude  of  the  signal  in  relative  units. 

The  absence  of  high  mode  in  frequency  spectrum  gives  an  opportunity  to  be  limited  by 
the  linear  vibrations  of  the  elastic  plate  only. 

By  the  use  of  SLM  the  following  value  of  the  decay  time  coefficient  at  ultrahigh  vacuum 
conditions  has  been  obtained 

p  =  (0.020  ±0.002)  sec _1. 

The  measurement  of  the  gas  pressure  in  the  vacuum  chamber  is  an  important  step  of 
the  measuring  technique.  Unfortunately,  construction  features  of  the  elastic  silicon 
single-crystal  plate  do  not  permit  to  realize  the  measurements  in  the  pressure  range 
upper  than  10'2  Pa.  The  reason  is  connected  with  a  small  damping  time  that  leads  to 
difficulties  in  the  decay  time  coefficient  determination. 

The  change  in  residual  gas  pressure  at  switched  off  means  of  gas  evacuation  provides 
natural  opportunity  to  carry  out  experiment  and  to  get  gas  pressure  dependence  of 
measuring  parameters. 
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During  the  total  time  of  experiment  of  more  than  three  months  the  gas  composition  of 
the  chamber  atmosphere  was  not  changed  significantly  and  is  consisted  of  96%  N2. 

At  the  gas  pressure  of  about  0.1  Pa  the  measurements  was  carried  out  with  the  rate  of 
1/15  musec-1  and  were  processed  by  FFT.  Frequency  spectrum  of  plate  vibrations 
excited  by  electrostatic  field  and  one  at  damping  conditions  are  presented  in  Fig. 3.6.  In 

this  case  the  decay  time  coefficient  gets  the  value  of  P  =  (8-6 ± °-7)sec  . 


f,  Hz 


(a)  (b) 

FIG.  3.6.  Frequency  spectrum  of  plate  vibrations  excited  by  electrostatic  field  (a)  and 
one  at  damping  conditions  (b).  This  is  the  Style  for  Figure  Captions.  Center  this  text  if  it 
doesn’t  run  for  more  than  one  line. 

The  analysis  of  the  results  gives  the  frequency  shift  of  about  17  Hz  that  is  corresponds 
to  1 ,4%.  The  frequency  increase  agrees  with  theoretical  predictions.  The  quantitative 
estimation  is  not  available  at  this  stage  because  of  the  problem  of  absolute  gas 
pressure  measurements. 

The  obtained  results  approve  the  initial  assumptions  and  demonstrate  agreement 
between  theory  and  experiment.  The  direction  of  frequency  shift  coincides  with  the 
theory  prediction. 

The  quantitative  estimation  problem  can  be  solved  by  geometric  parameters 
optimization  technique. 

The  problem  of  reception  of  quantitative  estimations  can  be  solved  by  change  of  the 
geometrical  sizes  of  an  elastic  plate.  This  decision  is  closely  connected  to  fabricating 
technique  of  Silicon  single-crystal  elastic  plate.  MEMS  technology  achievements  can  be 
useful  and  fruitful. 

The  other  way  of  decision  of  the  problem  is  connected  with  pressure  data  acquisition. 
Particularly  it  can  be  achieved  by  application  of  new  device  such  as  a  dynamic 
technique. 

Due  to  relative  simplicity  combined  with  high  efficiency  variation  methods  are  widely 
used  for  studying  an  elastic  body  deformation.  In  order  to  write  down  the  equation 
describing  process  of  the  plate  movement  in  a  gas  it  is  convenient  to  construct 
functional  of  the  moving  plate  and  determine  its  extreme  by  the  use  of  variation  Ritz 
method. 

Functional  variation  equation  written  for  whole  volume  of  the  elastic  body  has  the 
following  form 
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(3.5) 


8\\{p—  +  U (X))dvdt+  SB  =  0  , 

to  V  2 

SB  =  J  (p" )'  Swpla  +  J  pFiSWjdv  -  J psST dv , 


(3.6) 


where  -potential  energy  of  volume  element  of  the  body,  ^-material  density,  t() ,  1  - 

initial  and  present  time  respectively,  ^  -velocity  of  the  volume  element,  ^w,(X)_ 
displacement  deviations,  ^-volume  element,  $F -temperature  deviations  of  elasticity 

body,  F  -mass  force  applied  to  the  body,  5  -entropy  of  volume  element,  x  f 

Q  -  surface  area, 

(pn)i=pijnj, 

here  n>  -  normal  vector  of  J  axis,  Pij  -  stress  tensor. 

Taking  into  account  that  body  temperature  is  constant  last  item  in  expression  for  ^can 
be  discarded.  Some  principal  components  of  P,J  and  ^(X)  must  be  determined. 

In  this  study  the  free  damping  oscillations  of  the  square  thin  plate  with  rigidly  fixed 
edges  are  considered.  The  plate  is  placed  in  the  rarefied  gas.  Free  molecular  gas 
conditions  are  provided  by  high  ratio  of  mean  free  path  of  gas  molecules  to  linear  plate 
dimensions. 

Accepting  small  curvature  of  the  plate  for  the  first  degree  of  approximation  one  can 
neglect  tangential  stress  to  the  plate  surface. 

According  to  the  assumption  regarding  smallness  of  the  plate  curvature  the  normal 
vector  to  any  element  of  the  surface  can  be  considered  as  parallel  to  z-axis  (Figure  3.7). 
Inside  the  moving  plate  there  is  so-called  “neutral  layer”,  where  any  tensions  and 
compressions  are  absent.  The  neutral  layer  point’s  movement  is  described  by 
transverse  displacement  of  the  in-plate  co-ordinates  x  and  y  as  ® . 

Y 


FIG.  3.7.  Profile  the  neutral  layer  point’s  displacements  of  the  plate. 


li 


Therefore,  the  vector  of  displacements  of  points  of  the  neutral  layer  can  be  presented 
as 

co  =  (0,0,  c o0 ) 

5 

where  co°~  ^ x ’ y,t^ .  For  small  displacement  the  points  of  the  neutral  layer  move  in 
perpendicular  directions  relatively  to  the  plane  of  initial  plate  position. 


Let  us  introduce  the  coordinate  system  x  ,y  ,z  connected  with  the  moving  plate 
surface. 


If  to  take  into  account  that  unit  normal  to  any  surface  element  is  parallel  to  z-axes  the 
components  of  molecule  velocities  can  be  written  as 


< 


3  =3. 

3,=3.±*^ 
z  dt 


(3.7) 


8a>0 

The  different  sign  at  &  connected  with  the  opposite  plate  sides. 


The  normal  momentum  flux  to  the  elementary  surface  area  can  be  expressed  as 

Pi  =  (Bexp(-B2)  +  7ryH±-  +  B2)(l  +  erf(B )) , 


where  ; 


B  =  ± 


(3.8) 


1 


da>n 


-mP,  8<  p. 


and  -  density  and  temperature  of 


undisturbed  gas  correspondingly,  R  -  gas  constant. 


According  to  accepted  Knudsen  definition  of  accommodation  coefficient  the  equation  for 
normal  momentum  flux  created  by  reflected  molecules  can  be  presented  as 

Pr=0-~an)Pi+anPw  (3'9) 


The  expression  for  total  normal  momentum  flux  applied  to  the  plate  has  a  view  as 

P  =  Pi+Pr=(2~ocn)pi+anpw  (3.10) 

where  GCn  -  so-called  normal  momentum  accommodation  coefficient,  Pi  -  normal 

momentum  flux  created  by  incident  molecules,  Pw  -  normal  momentum  flux  caused  by 
the  influence  of  reflected  molecule  having  Maxwell  distribution  function  with  zero 
macroscopic  velocity  and  with  the  temperature  of  the  surface. 

Assuming  velocity  loss  of  diffuse  reflected  molecules  that  leave  the  surface  with 
equilibrium  distribution  function  one  could  write  Pw  as 

p.=^2nBT„N„  (3.11) 

where  m  ■  molecule  mass,  '  -  number  of  molecules  striking  the  unit  surface  per  unit 
time: 
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(3.12) 


N,  =  n 


I  co  i 


\KT\ 

In 


(exp(-5  )  +  ~JnB(\  +  erf(B)) , 


where  -  gas  numerical  density  at  infinity. 

Oscillations  are  believed  to  be  at  quasi-steady  condition.  Hence,  it  can  be  considered 
that  B2  « 1. 


Finally,  the  normal  pressure  to  the  surface  unit  becomes 

P  =  Px±PxJ^^(anj  +  2-an), 

V  n  ot  4 

p<°  -  gas  pressure  at  infinity,  P«  =  nPn . 

Elastic  body  presents  the  plate  of  cubic  syngony  with  (100)  orientation. 


(3.13) 


Density  of  strain  energy  u  in  Hookian  approximation  is  a  quadratic  function  of 
components  of  the  deformation  tensor: 


1  6  6 
Z  A=1  ju= 1 


A? ju  ’ 


(3.14) 


where  e*,  fl  -  deformation  tensors,  ^ ,  P  indexes  are  determined  as 

1  =  xx  ,  2  =  yy  ,  3  =  zz,  4  =  yz ,  5  =  zx ,  6  =  xy  , 

^  -  components  of  the  elasticity  tensor  (cubic  syngony)  can  be  written  in  a  matrix 
form: 


r. 


C,  = 

¥ 


'll 


'12 


11  12 

C, 


0  0 


•a 

0 


0 

0 

0 

C44 


0  0  0  0 
0  0  0  0 


0 

0 

0 

0 

C44 


0  ^ 

0 

0 

0 

0 


0  c 


44/ 


Hence, 


U  =  —  c,  ,(e2  +  e2  +  e°  )  H — ciA(e 7  +  e2  +  e7  )  +  c,.,(e  e  -fee  -fee).  (3.1 5) 

2  11''  yy  ZZ '  ^  44  ^  yz  xy-'  12V  yy  zz  zz  x*  x*  ;yy  ^  \  / 


Taking  into  account  assumptions  made  above  and  using  correspondent  expressions 
one  can  obtain 


erv  =  -z- 


d'COr. 


d2a> n  .  d2con 

W'e*’  z^'ea=e*=e* 


(3.16) 


Expression  for  U  can  be  obtained  by  replacing  (3.16)  in  (3.15) 

1  „  2,32<ynx2  d2co0  d2co0 


(3.17) 


Finally  for  the  whole  plate  it  becomes: 
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av,„  )2 + ^  +  2Cjj(a^)2 + cJpydyy)dxdy  _ 


5x2 


dy 


dxdy 


dx 1  dy 1 


(3.18) 


Present  consideration  of  the  variation  task  is  based  on  following  assumptions. 

a)  The  damping  of  the  plate  oscillations  caused  by  internal  stresses  is  negligible  in 
comparison  with  the  damping  due  to  external  forces  due  to  gas-surface 
interactions. 

b)  The  plate  produces  the  small  oscillations  in  strongly  perpendicular  directions  to 
the  plane  of  unmoved  plate. 

c)  The  oscillations  are  believed  to  have  a  single  harmonic. 

Equations  (2.5)  and  (2.6)  are  transformed  to  the  following  form: 

5 J (T  +  U)dt  +  2j  jj  px  J2^7-  ^  (an  -  +  2  -  an  )5coadxdydt  =  0 ,  (3.19) 

<0  h  n  *  n  dt  4 


where  T -kinetic  energy  of  the  plate  determined  as 

T  =  ^\\  i^-fdxdy,  (3.20) 

2  "  dt 

here  h  -  plate  thickness. 

Determination  of  the  extreme  of  the  functional  (2.19)  is  made  by  the  Ritz  method. 

The  minimization  form  is  taken  in  the  following  view: 

coQ(x,y,t )  =  2^anm(t) sin  - sin-  - —  ,  (3.21 ) 

n,m  d  d 

where  a  -  side  length  of  the  plate,  n  ,m  e 

This  form  must  satisfy  to  boundary  conditions  corresponding  to  the  case  of  the  plate 
with  rigidly  fixed  edges. 

coQ  =  0  at  x  =  0,  x  =  a  and  y  =  0,  y  =  a ; 

=  0  at  x  =  0,  x  =  a  ;  ^SL  =  o  at  y  =  0,y  =  a  . 
dx  dy 


In  addition  to  the  conditions  mentioned  above,  function 


a>0(x,y,t) 


in  continuous  region 


G  e  [0, a] x  [0, a]  at  n  1  ,m  1  ^as  a  sjng|e  extreme  point  in  the  center  of  the  plate 
which  corresponds  to  the  oscillation  condition  with  single  harmonic  and  the  basis 


a  a 
2’ 2 


.  2  Hxn  ■  2  nym 
sin  srn-^ — 

functions  a  and  a  have  the  spanning  property. 


Using  correspondent  replacements  and  some  lemmas  of  variation  calculus,  one  can 
obtain  the  definitive  equation 

da(t)  32  /?  V 

- 1 - —  (3cn  +2c22  +  cl2)a(t)  =  0  .  (3.22) 


,  2  d2a(t) 

Pha  ,2  +  2P« 


n 


dtA 


2 RT ,  „  „ 

- («  -  +  2-a„)fl- 

n  4 


dt 


21a 


Finally,  the  damping  coefficient  of  plate  oscillations  can  be  presented  as  a  function  of 
normal  momentum  accommodation  coefficient  by  the  following  way 
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(3.23) 


px  2RT  n 

P  =  ^fJ - («,t  +  2  ~a„). 

ph  V  n  4 


This  expression  can  be  compared  with  the  result  obtained  earlier: 

P  =  ^^(2~an). 

ph  V  n 


(3.24) 


The  normal  momentum  accommodation  coefficient  dependence  of  P  in  relative  units 
obtained  from  equations  (3.23)  and  (3.24)  is  shown  in  Figure  3.8. 


On 


Fig.  3.8.  Normal  momentum  accommodation  coefficient  dependence  of  P  related  to 

P  / 

/ Ph  for  equations  (3.23)  and  (3.24). 


The  contribution  of  normal  momentum  accommodation  in  damping  coefficient  according 

to  equation  (3.23)  is  about  20%  at  The  difference  between  P  obtained  from 
equations  (3.23)  and  (3.24)  can  reach  up  to  78%. 


Thus,  equation  (3.23)  is  more  preferable  to  processing  experimental  data  especially  for 
the  gases  with  normal  momentum  accommodation  coefficient  close  to  one. 


YEAR  2: 


TASK  4:  ENERGY  ACCOMMODATION  MEASUREMENTS. 


Figures  4.1  a  and  4.2b  show  the  measured  power  supplied  to  the  tungsten  sample  in 
Argon  and  CO2,  respectively.  Each  symbol  represents  the  average  of  more  than  100 
data  points  acquired  by  the  computer  controlled  heating  and  temperature  stabilization 
system  operating  at  an  input  sampling  rate  of  50  kHz. 
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q/q\ 


Fig.  4.1 .  Electric  power  supplied  to  sample  vs.  temperature 
for  (a)  argon,  (b)  C02,  and  (c)  at  discharge  current  J  = 

(1)0,  (2)  40,  (3)  1 00,  and  (4)  200  mA. 

To  facilitate  analysis,  the  data  are  represented  as  the  ratios  q/qi  between  the  powers 
measured  at  particular  temperature  and  T  =  1000  K.  The  curves  obtained  for  different 
gases  exhibit  different  behavior,  particularly  at  high  temperatures. 

The  measurement  parameter  settings  are  such  that  the  relative  error  does  not  exceed 
5%. 


The  power  supplied  to  the  glow  discharge  also  has  an  absolute  error  due  to  the  filament 
sputtering  during  the  measurements  performed  at  operating  values  of  current  and 
voltage.  This  is  supported  by  the  results  of  preliminary  scanning  electron  microscopy 
studies.  Additional  indirect  evidence  is  provided  by  the  fact  that  the  sample  resistance 
increases  by  approximately  10%  over  a  measurement  cycle,  which  corresponds  to  a 
decrease  in  wire  radius  by  up  to  3%. 

We  have  also  performed  experiments  on  glow  discharge  in  air  (see  Fig. 4. 1c).  Their 
results,  being  similar  to  those  obtained  for  C02  plasma,  suggest  that  the  characteristics 
of  heat  transfer  at  high  temperatures  strongly  depend  on  the  composition  of  the 
molecular  gas  plasma  produced  in  the  glow  discharge.  The  analysis  that  follows 
focuses  on  heat  transfer  in  Argon  and  C02  plasmas. 


The  most  important  plasma  processes  contributing  to  heat  transfer  across  the  tungsten 
surface  include  radiation,  sputtering,  molecular  heat  conduction,  and  heterogeneous 
recombination. 
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The  radiative  heat  fluxes  are  approximately  equal  for  both  gases.  Indeed,  even  if  a  layer 
of  some  compound  forms  on  the  filament  surface  in  C02  atmosphere,  the  ensuing 
change  in  emissivity  can  be  neglected  in  view  of  the  concurrent  sputtering;  i.e.,  the 
radiative  heat  flux  can  be  assumed  to  depend  only  on  the  filament  temperature. 

Neither  physical  sputtering,  nor  chemical  sputtering,  nor  sublimation  can  be  singled  out 
as  a  dominant  process  responsible  for  material  loss  when  comparing  the  two  gases 
under  study.  Furthermore,  the  sputtering  rate  constant  is  independent  of  sample 
temperature  far  from  the  sublimation  point,  and  sputtering  has  similar  characteristics  in 
argon  and  C02  atmospheres  when  all  other  parameters  are  equal. 

Since  the  Knudsen  number  for  the  actual  geometry  used  in  our  experiments  is  Kn  >  1 0, 
analysis  can  be  restricted  to  free  molecular  heat  transfer.  The  gases  under  study  have 
similar  molecular  masses,  and  preliminary  estimates  give  a  degree  of  ionization  varying 
between  10“5  and  10“3.  Therefore,  the  thermal  accommodation  coefficients  for  Ar-W 
and  C02-W  systems  are  approximately  equal  according  to  Baule’s  theory. 

Based  on  the  discussion  above,  we  suggest  that  the  experimentally  found  difference  in 
heat  transfer  between  the  gases  can  be  attributed  to  heterogeneous  recombination 
processes,  which  are  significant  at  a  temperature  of  1700  K. 

Consider  the  normalized  supplied  power  q/ql  plotted  as  a  function  of  discharge  current 
for  C02  and  Argon  in  Fig. 4. 2.  The  slopes  of  the  lines,  obtained  by  least  squares 
interpolation,  differ  by  a  factor  of  more  than  2. 


<7/<7l 


Fig.  4.2.  Electric  power  supplied  to  tungsten  wire  vs.  dis¬ 
charge  current  at  sample  temperature  T  =  1700  K:  (1) 
argon;  (2)  C02. 

In  the  case  of  argon,  the  linear  interpolation  slopeis  determined  by  the  heating  of  the 
discharge  chamber  walls  throughout  the  measurement  process,  which  can  be  caused 
both  by  radiation  from  the  tungsten  wire  and  by  heat  transfer  from  the  plasma.  These 
effects  can  be  taken  into  account  by  introducing  certain  changes  into  the  measurement 
procedure. 

Glow  discharge  in  carbon  dioxide  is  characterized  by  complex  plasma  composition 
strongly  depending  on  discharge  current.  The  dominant  mechanism 
of  dissociation  is  the  following: 

CO2  +  €  — ►  CO2  — *■  CO  +  o. 


Its  products  can  obviously  play  an  important  role  in  plasma-surface  interaction. 
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Heterogeneous  recombination,  which  significantly  contributes  to  the  heat  transfer 
between  the  tungsten  filament  and  the  C02  plasma,  may  follow  two  mechanisms:  the 
Eley-Rideal  mechanism  (reaction  between  a  gas  phase  species  and  an  adsorbed  one) 
and  the  Langmuir-Hinsheiwood  mechanism  (reaction  between  adsorbed  molecules). 

A  general  analysis  may  include  a  wide  variety  of  recombination-dissociation  and 
exchange  reactions.  However,  the  available  experimental  data  are  not  sufficient  to 
single  out  a  dominant  mechanism. 

Our  experimental  results  demonstrate  that  the  hot  filament  method  can  be  used  under 
glow  discharge  conditions  to  examine  characteristics  of  heat  transfer  between  low 
temperature  monatomic  or  molecular  gas  plasma  and  a  solid  surface. 

A  qualitative  analysis  of  the  measured  results  suggests  that  heterogeneous 
recombination  plays  a  significant  role  in  plasmas  of  C02  and  other  polyatomic  gases, 
such  as  those  contained  in  air. 

Task  5:  Surface  structure  effects  on  gas  flow  rate  in  channels. 

One  of  the  factors  that  have  a  great  influence  on  gas/surface  interaction  and  as  a  result 
on  heat  and  mass  exchange  in  a  gas/surface  system  is  surface  roughness.  This  effect 
has  been  studied  experimentally  earlier  and  described  with  the  use  of  different  surface 
structure  models  later.  The  main  problem  of  surface  structure  modeling  at  least  at  the 
first  steps  was  a  lack  of  the  information  on  features  of  real  surfaces.  Significant  progress 
in  experimental  technique,  in  a  particular,  numerous  realizations  of  Scanning  Probe 
Microscopy  open  new  opportunities  in  surface  diagnostics  at  micro-  and  nano-  scales. 
Because  of  essential  increase  in  a  role  that  gas/surface  interaction  plays  at 
miniaturization  of  gas  dynamics  systems  the  problem  becomes  rather  topical. 

In  this  study  an  attempt  to  describe  gas  particles  scattering  by  the  real  Si-based  surface 
has  been  made  with  the  use  of  DSCM  approach.  The  surface  morphology  of  the  sample 
has  been  studied  by  Atomic  Force  Microscopy  (AFM)  that  provides  obtaining  such 
characteristics  of  the  surface  as  mean  square  height  of  micro  roughness,  fractal 
dimension  and  micro  roughness  distribution  function.  The  surface  of  the  sample  studied 
by  AFM  presents  area  of  20x20  mkm  and  characterizes  by  160  thousand  points.  The 
surface  has  been  simulated  by  triangulations  which  are  necessary  to  calculate 
coordinates  of  crossing  points  of  particle  trajectories  with  a  surface.  Fragment  of  the 
surface  simulated  according  to  AFM  data  is  presented  in  Fig. 5.1 . 

DSMC  in  a  probe  particle  version  has  been  realized  for  modeling  of  gas  scattering  by  a 
surface  at  various  parameters  of  micro  roughness  considering  diffuse  or  specula 
reflection  of  gas  particles  from  a  separate  elementary  plate  of  a  rough  surface.  As  a  first 
step  monochromatic  gas  particles  flow  directed  to  the  surface  at  various  angles  is 
considered.  The  next  step  supposes  to  consider  gas  flow  in  a  channel  and  take  into 
account  real  distribution  function  of  micro  roughness  and  gas  molecule  velocity 
distribution  function  on  an  input  in  the  channel  as  well. 

As  a  result  of  the  study  index  ellipsoids  (scattering  indicatrix)  for  various  angles  of 
incident  particles  and  parameters  of  surface  roughness  have  been  obtained.  Lobular 
and  completely  diffuse  scattering  presented  in  Fig. 5. 2  demonstrate  effect  of  variation  in 
surface  roughness  factors.  Test  calculations  confirm  well  known  fact  of  completely 
diffuse  scattering  of  monochromatic  particles  flow  with  thermal  energy  by  a  rough 
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surface  and  specula  reflection  in  the  case  of  a  smooth  surface.  It  has  been  found  that 
number  of  “particle  -  separate  rough  surface  side”  collisions  depends  on  surface 
roughness  parameters  and  as  a  result  influences  on  scattering  indicatrix.  Fig. 5. 3 
illustrates  velocity  distribution  function  for  reflected  particles  simulated  on  the  basis  of 
100  thousand  tests. 


Fig. 5.1.  Surface  structure 
simulated  by  triangulations 
according  to  AFM  data 


Fig. 5. 2.  Scattering  indicatrix 
as  a  function  of  average 
relative  height  of  surface 
roughness:  ■,  .,*-1,2, 
and  >  3,  respectively 


Fig. 5. 3.  Gas  velocity 
distribution  function  for 
reflected  particles  (y- 
component  in  relation  to 
most  probable  value) 


One  of  questions  of  principle  which  was  put  at  this  stage  consisted  in  finding  out  how 
much  gas  flow  rate  is  sensitive  to  character  of  surface  structure  at  the  same  size  of 
average  height  of  microroughness.  The  thing  is  in  some  cases  the  channels  have 
specially  prepared  surface  with  regular  structure.  In  other  cases  we  have  real  surfaces 
with  statistical  distribution  of  microroughness. 
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■  -  -  Cl au sing's  data 


- Test  calcdation 


— Channel 
according  to 
AFM 

“Statistical” 

surface 


— •—  Regular  surface 
structure 


01  23456789  10 

UR  (channel  lenc^it  to  radius  ratio) 

Fig. 5. 4.  Calculated  results  for  passing  probability  of  cylindrical  channels  with  various 

surface  structures. 

As  it  follows  from  Fig. 5. 4  there  is  an  essential  distinction  between  the  results  obtained 
for  smooth  surface  (Clausing’s  data)  and  for  surfaces  with  regular  or  statistical  structure 
as  well  as  for  surface  reconstructed  according  to  AFM  data.  It  means  that  the  use  of 
simplified  model  of  a  surface  as  it  realized  in  many  cases  can  bring  considerable  errors 
in  calculation  of  gas  flow  rate. 

Task  6:  Surface  chemical  composition  effects  on  molecular  heat  transfer. 

The  experimental  setup  created  for  realization  of  this  stage  provided  measurements 
with  ultrahigh  vacuum  conditions  up  to  1 0"10  mmHg  as  well  as  with  means  of  surface 
and  gas  phase  diagnostics  by  the  use  of  vacuum  system  based  on  the  helium  cryogen 
pump,  mass-spectrometer  and  Auger-electron  spectrometer  built  in  measuring 
chamber. 

The  situation  observed  in  measuring  system  during  surface  purification  process  and 
thermal  accommodation  measurements  in  He-W system  is  presented  in  Figs. 6. 1-6.3. 


It  units 


2  4  H 16 18  28  40 

Atomic  mass 


Fig.  6.1 .  Mass  spectrum  in  vacuum  system  before  (left)  and  after  (right)  leak-in  He. 
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Fig. 6. 2.  Energy  peaks  of  Auger-electrons  at  the  process  of  tungsten  surface  covering  by 
hydrocarbons:  bottom  peaks  -  “clean”  surface;  top  peaks  -  “dirty”  surface. 


dN/dE 


Fig. 6. 3.  Auger  spectrum  for  W wire  sample  at  different  stages  of  cleaning:  1  -  surface 
after  series  of  purification  cycles;  2-  partially  contaminated  surface  with  C,  N,  Oatom 

complexes. 

The  intensity  of  the  peaks  in  Fig. 6. 2  and  Fig. 6. 3  corresponds  to  number  of  Auger- 
electrons  that  reflects  concentration  of  atoms  on  the  surface.  The  term  “dirty”  surface 
means  here  the  surface  at  initial  conditions  before  applying  special  purification 
measures. 

The  value  of  thermal  accommodation  coefficient  aE  at  the  process  of  surface 
contamination  removal  changes  from  ~  0.45  for  the  surface  without  special  treatment 
(typical  value  for  “dirty”  surfaces)  to  -0.02.  Purification  cycles  included  annealing  at 
ultrahigh  vacuum  conditions  and  chemical  reduction  in  hydrogen  environment. 
Dependence  of  aE  on  surface  temperature  and  surface  coverage  9  (surface  fraction 
covered  by  adsorbed  gases)  for  He- W system  is  presented  in  Fig. 6. 4. 
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Fig. 6. 4.  Thermal  accommodation  coefficient  cte  as  a  function  of  6  and  surface 

temperature: 

7  -298  K\  2-446  K\  3-968  K\  4-  1650  K\  5-2020  K. 

On  the  basis  of  classical  knowledge  about  movement  of  atoms  and  lattice  theory  of  F. 
Goodman  and  G.  Wachman  the  program  modeling  equilibrium  and  nonequilibrium 
scattering  of  helium  atoms  on  a  three-dimensional  crystal  tungsten  lattice  taking  into 
account  partial  surface  coverage  with  adsorbed  oxygen  atoms  is  developed.  The 
modeling  scheme  is  presented  in  Fig. 6. 5. 


2*  /•> 


Fig. 6. 5.  Scattering  of  He  atom  on  W  lattice  partially  covered  by  oxygen  atoms 
due  to  adsorption,  6  ~  0.8:  1,2-  positions  of  He  atom  (initial  and  finite,  respectively). 

The  scattering  model  is  constructed  on  the  basis  of  following  assumptions:  a)  atoms  of 
gas  and  a  solid  body  represent  spheres,  their  radius  and  mass  correspond  to 
characteristics  of  real  atoms;  b)  W  sample  is  modeled  as  crystal  with  bcc  lattice 
structure  producing  thermal  vibrations  at  temperature  Ts;  c)  forces  of  interaction 
between  atoms  of  a  crystal  are  modeled  only  as  ones  between  the  nearest  neighbors; 

d)  particles  of  gas  falling  on  a  surface  represent  the  Maxwell’s  flow  at  temperature  Tg; 

e)  interaction  force  between  atom  Fie  and  a  solid  body  is  calculated  taking  into  account 
all  atoms  of  a  crystal  and  the  atoms  adsorbed  on  it;  f)  adsorbate  is  modeled  by  atoms 
of  oxygen  as  their  primary  presence  on  a  surface  proves  to  be  true  by  the  data  of  the 
electron-spectroscopic  analysis;  g)  the  case  of  thermal  scattering  of  atoms  on  a  crystal 

surface  with  <  0.1eV  is  considered  as  far  as  scattering  of  particles  with  higher 
energy  demands  quantum-mechanical  interpretation;  h)  influence  of  every  possible 
directions  of  crystal  axes  on  scattering  process  is  not  considered,  otherwise,  atom  of 
gas  interacts  with  a  flat  surface.  Thus,  the  data  obtained  during  modeling  represents 
integrated  behavior  of  gas  -  surface  system  that  is  realized  during  all  possible 
processes  of  scattering. 
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Interaction  of  helium  atom  with  a  solid  body  is  described  by  potential  function  of  pair 
interaction  between  gas  atom  and  each  atom  of  a  solid  body  (including  the  adsorbed 
atoms).  As  interaction  potential  between  atoms  of  all  types  the  Morze  potential  has 
been  chosen. 

In  terms  of  energy  accommodation  coefficient  aE  an  efficiency  of  molecular  heat 
exchange  of  helium  with  clean  and  partially  covered  by  adsorbed  oxygen  atoms 
tungsten  surface  is  calculated  for  various  temperatures.  The  calculated  temperature 
dependences  are  compared  with  the  data  obtained  in  experiments  with  surface 
diagnostics.  Within  the  limits  of  the  developed  approach  the  results  of  calculations  are  in 
a  satisfactory  agreement  with  experimental  data  obtained  for  equilibrium  aE  and  clean 
surface  (Fig. 6. 6.)  as  well  as  for  nonequilibrium  aE  and  a  surface  partially  covered  by 
oxygen  atoms  (Fig. 6. 7). 


Fig. 6. 6.  “Equilibrium”  aE  for  clean  surface. 
Simulation  (bold  line)  and  experiment: 

1  -  L.B.  Thomas,  1967 
2-  J.  Kouptsidis  et  al,  1970. 


Fig. 6. 7.  ’’Nonequillibrium”  arEfor 
contaminated  surface  Tg=  300K\ 
Ts=298K\  1-  simulation, 

2-  experiment, 
Ts=446K',  3-  simulation, 

4-  experiment. 


Task  7:  Surface  chemical  composition  effects  on  freemolecular  gas  flow  in 
channel. 

In  our  case  the  surface  preparation  procedure  is  based  on  evaporation  and  deposition 
of  metal  on  internal  surface  of  a  cylindrical  channel.  As  the  investigated  channel  the 
Pyrex  glass  cylindrical  channel  in  diameter  of  6.4  mm  and  length  150  mm  has  been 
taken.  Titanium  filament  with  0.6  mm  in  diameter  is  stretched  along  its  axis.  Heating  of 
the  wire  is  reached  by  the  use  of  electric  current.  The  change  in  chemical  composition 
of  the  surface  is  controlled  by  Auger-electron  spectroscopy  using  special  pilot  sample 
(witness)  that  is  situated  near  the  channel  and  simultaneously  exposed.  To  arrange  the 
guaranteed  surface  coverage  with  adsorbate  a  special  oxygen  source  has  been  used. 

The  results  of  surface  analysis  for  two  situations  (“dirty”  and  “clean”  surface)  are 
presented  in  Fig. 7.1 .  Similar  spectrum  in  vacuum  system  one  can  observe  practically  in 
all  cases  when  special  measures  for  purification  are  not  applied.  This  result  is  expected 
if  take  into  account  mass-spectrum  of  typical  vacuum  system  (Fig. 6.1 ,  left). 


23 


E,  eV 


Fig. 7.1 .  Auger-spectrum  for  “clean”  (top)  and  “dirty”  (bottom)  surface. 

Actually  the  real  surface  is  oxide  covered  by  adsorbed  gases  from  contaminations  of 
investigated  gas  (for  instance,  helium)  and  residual  environment  of  vacuum  system. 
Both  adsorption  and  chemisorptions  take  place.  As  we  see  at  the  lower  spectrum 
(Fig. 7.1 )  there  are  some  picks  belonged  to  carbon,  nitrogen  and  oxygen.  It  is  a  typical 
situation  as  far  as  any  vacuum  system  has  primary  hydro  carbonic  complexes,  oxygen, 
nitrogen  and  its  complexes.  All  of  them  participate  in  formation  of  surface  covering. 
Hydrogen  we  cannot  see  by  Auger-electron  spectroscopy  but  it  is  present  in 
hydrocarbons  that  are  identified  as  vacuum  grease.  It  explains  why  such  surface  is 
considered  as  “dirty”. 

An  experience  shows  that  at  the  first  stage  of  purification  procedure  some 
contaminations  can  be  removed  rather  easily  by  heating  of  the  sample.  But  it  not 
concerns  oxygen.  That  is  why  in  our  simulations  we  take  oxygen  atoms  as  main 
component  of  surface  coverage  that  still  is  present  on  the  surface  even  after  some 
purification  cycles. 

Other  way  to  have  guaranteed  surface  covering  is  to  fill  system  with  pure  oxygen  when 
the  surface  has  already  completely  cleaned  by  applying  special  surface  purification 
procedure.  As  it  follows  from  top  spectrum  in  Fig. 7.1  there  is  only  Auger-pick  from 
titanium  and  nothing  else.  It’s  clean  surface.  Though  hydrogen  we  can  not  see,  it  is  not 
present  there  as  it  leaves  the  surface  after  heating  at  about  600  K. 

The  comparison  of  Auger-spectrum  obtained  for  tungsten  wire  and  titanium  vacuum- 
deposited  surface  shows  that  in  both  cases  the  surface  contamination  is  consisted  of 
the  practically  the  same  atom  complexes. 

The  results  of  free  molecular  gas  flow  rate  measurements  in  channel  as  a  function  of 
oxygen  exposition  is  presented  in  Fig. 7. 2. 
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Fig. 7. 2.  Mass  flow  rate  M'  in  a  cylindrical  channel  for  some  noble  gases  in  free 
molecular  regime  as  a  function  of  oxygen  exposition  (in  Langmuir).  Bold  line  -  Test 
Particle  Monte  Carlo  method  (TPMC-simulation). 

On  the  bases  of  Test  Particle  method  the  dependence  of  gas  flow  rate  on  titanium 
surface  coverage  with  oxygen  has  been  obtained.  An  original  surface  stricture  is 
simulated  by  hep  crystal  lattice  (unlike  tungsten  surface  that  is  simulated  by  bcc  lattice). 
The  process  of  gas  atom  scattering  by  the  surface  is  simulated  similar  to  the  procedure 
realized  in  the  case  of  thermal  accommodation  for  He- W system.  The  preliminary 
results  of  simulation  obtained  for  free  molecular  flow  rate  of  He  in  the  channel  with  clean 
and  completely  covered  by  oxygen  77  surface  are  close  to  experimental  results. 
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